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Abstract 

We investigate the quenching process in lattice systems with short range 
interaction and several crystalline states as ground states. We consider in 
particular the following systems on square lattice: 

- hard particle (exclusion) model 

- q states planar Potts model. 

The system is initially in a homogeneous disordered phase and relaxes toward 
a new equilibrium state as soon as the temperature is rapidly lowered. The 
time evolution can be described numerically by a stochastic process such as 
the Metropolis algorithm. The number of pure, equivalent, ground states is 
q for the Potts model and r for the hard particle model, and it is known that 
for r, or q > d+1, the final equilibrium state may be polycrystalline, i.e. not 
made of a uniform phase. We find that in addition Ug and qg exist such that 
for r > Tg OT q > qg the system evolves toward a glassy state, i.e. a state 
in which the ratio of the interaction energy among the different crystalline 
phases to the total energy of the system never vanishes; moreover we find 
indications that Vg = qg. We infer that q = qg (and r = Vg) corresponds 
to the crossing from second order to discontinuous transition in the phase 
diagram of the system. 
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Introduction 



A large variety of materials occurs as glass or amorphous (Elliot 1990) but 
mechanisms driving matter into such kind of states are far from clear. Be- 
cause of the inherent complexity of the problem, much effort has been de- 
voted also to the investigation of simple models hke lattice models. These 
systems may reproduce at least a few of the characteristic features of glassy 
behaviour, however they arc generally characterized by quenched disorder or 
competing interactions. In contrast with real glass former such model sys- 
tems are not able to crystallise and therefore cannot give any information on 
possible transitions from glassy to crystalline behaviour. 

In this work we show that also simple lattice models with short range 
interaction and no quenched disorder may enter a kind of glassy phase. In 
particolar we show in section 2 that if a g state Potts system is quenched 
from infinite to zero temperature, it relaxes toward a state of minimum energy 
only if q is small enough. Otherwise it slowly equilibrates into a state with 
energy well above the ground state. Analogous observations on different 
lattice models are elaborated in section 3 and indicate that such behaviour 
may be universal. 

The Potts model 

Consider a regular lattice in which each site can stay in one of q distinct 
states, or colours. To each site one may associate a variable rji that takes 
the values 0, 1, 2, . . . , g — 1. The energy of a configuration is given by (Potts 
1952) 

E^-e^{S,,,,,-l), (1) 

where the summation is over pairs of nearest neighbour sites. If two nearest 
neighbour sites have the same colour their contribution to the total energy 
vanishes, otherwise they contribute with an energy e. There are therefore q 
distinct ground states; in each of them all sites in the lattice have the same 
colour and the total energy takes the absolute minimum value, E = 0. These 
states are crystalline states. 

Let Up the number of site pairs such that the two sites have distinct 
colours, that is 

np = J2a-S,,,,^), (2) 

(ij) 
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thus 

E = e Up (3) 

and in a crystalline state = 0. A polycrystalline state is composed of 
clusters of crystalline states of distinct colours separated by Up site pairs of 
different colours. The number Up however must be small. More precisely we 
call a state polycrystalline when 1/ — Up/N — > as the number of lattice sites 
iV ^ 00, so that the total energy per site u — E/N — > 0, the same of the 
crystalline state. 

We have simulated a quench to zero temperature as follows (Metropolis 
et al. 1953). At each step of simulation a site is chosen at random and its 

new colour is chosen at random. If the energy does not increase, the new 
colour is taken as the new state. The initial configuration of the system is a 
high temperature one, that is each site has a colour chosen at random. 

We have studied the time behaviour of the energy per site u = E/N 
after the quench on a square lattice for several values oi N — L'^. It is well 
known that the dynamics of the Potts model usually gets pinned very soon 
into a non equilibrium state because of the finite size of lattices used in the 
simulations (Glazier et al. 1990). In order to escape this trap we have used 
both periodic boundary conditions and fixed boundary conditions. In the 
latter, all sites on the lattice boundary are frozen in a given state, say the 
r) = state. By this way one can discriminate finite size and boundary effects 
in the behaviour of the system. In fact we have assumed that, when plotting 
energy vs time, the siipcrimposing section of curves resulting from different 
simulation conditions, may be identified with the genuine relaxation of the 
system. We have extrapolated the behaviour for L — > 00 from the l/^JJt) 
part of these curves. As examples, the plots resulting from q — 3 and q — 7 
are shown in figures 1 and 2. 

From simulations on a square lattice we have drawn the following conclu- 
sions: 

a) for g = 2 the system relaxes toward a crystalline state: Up ^ and 

li~rV2; (4) 

b) ioT q — 3 the state is polycrystaUine: Up ^ but u — Up/N — > 0, 

u ~ r^^. (5) 
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This agrees with the Lifshitz proposition (Lifshitz 1962) that the final states 
are non pure crystalhne when q > d + 1 {q > 3 in d = 2). However the 
proposition does not say anything about the energy of these states; 

c) for > 4 we observe that u ^ u* > when t ^ oo, since u ^ u* > 0, i.e. 

The final state is not polycrystalline and, in addition, its energy is above the 
energy of the ground state (in fact it is seen to be a monotonic function of 
q). We call such pinned state a glassy state. 

d) for g = 4 the asymptotic behaviour is well fitted by 

u^t-^/^\n{t), (7) 

Preliminary results from simulations of quenches at T > indicate that the 
observed features persist as far as T is not too large. 

Hard core particles on a lattice 

Extension of the above observations to difi^erent lattice models, also char- 
acterized by short range interaction and absence of quenched disorder, leads 
to remarkable conclusions. We discuss here some results reported in the liter- 
ature for a class of lattice gas models, namely exclusion models with random 
sequential adsorption and in-plane diffusional dynamics. 

These models consist of a regular lattice in which particles are placed at 
the sites of the lattice. The deposition of a particle in a site blocks that site 
and the sites in a certain neighbourhood, whereas the energy is simply given 

by 

E = -ii{n - Ucp) (8) 

where > is the chemical potential, n is the total number of deposited 

particles and n^p is the number of deposited particles in the closest packed 
configuration. The configuration of closest packing is a ground state config- 
uration because the number of particles gets maximum so that the energy 
E is minimum (/x > 0) (in analogy with the Potts model the ground state 
energy is arbitrarily set to zero). 

If a particle blocks the deposition of other particles in a certain neighbour- 
hood the ground state is not unique. When for instance nearest neighbours 
sites are blocked on a square lattice (NN model) then there are two ground 
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states corresponding to place the particles on either one or the other of two 
sub-lattices. If a particle blocks the nearest neighbours as well the next near- 
est neighbours (NNN model) then there will be four relevant ground states, 
corresponding to four different sub-lattices. If besides nearest neighbours 
and next nearest neighbours the particle blocks also the next next nearest 
neighbours (NNNN model) then the number of ground states will be five. 

Once deposited, particles may diffuse along the lattice, with the prescrip- 
tion that they cannot overlap other particles. This process, called diffusional 
relaxation, makes possible to depose more particles and, in principle, the 
system enter an ordered phase, coming eventually to the closest packed state 
where density cannot further increase. 

Depending on the model (NN, NNN, NNNN, etc.), the final state may be 
crystalline, polycrystalline or glassy, and can be classified in a way similar to 
the case of the Potts model discussed above. Let be the number of sites 
at the surface of clusters of sites in the same state, i.e. occupying the same 
sub- lattice: 

Us = ricp - n. (9) 

In analogy with the Potts model, if = the state is monocrystalline. If 
Us ^ but Ug/N ^ as iV ^ oo, the state is polycrystalline, and is glassy 
in the opposite case. 

We briefly elaborate here on some results reported in the literature by 
Wang, Nielaba and Privman (1993, 1993a), and by Eisenberg and Baram 
(1998, 2000). In the simulation a lattice initially empty is filled sequentially 
by placing particles at the non excluded sites. Deposited particles may diffuse 
to next neighbour sites if these are not blocked or occupied. Although the 
implementation of this dynamics varies slightly from one author to the other, 
their results are of general validity. Let us denote by r the number of ground 
states so that ricp = N/r, and let u = E/N. By analysing the papers cited 
above one concludes that: 

a) for r=2 (NN model) — > and m — >^ (Wang et al. 1993) since 

Us ^ (iQ) 

b) for r = 4 (NNN model), Us ^ Q but Us/N ^ as ^ oo, (Wang et 
al. 1993a, Eisenberg and Baram 1998) and w ^ as 

Us ^ ln(t)r^/2. (11) 
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c) for r = 5, M — > > and 

Us^r^ + ul; <>0. (12) 

This behaviour has been observed in the very late stage (Eisenberg and 
Baram 2000), when the (finite) system is almost completely filled, and there- 
fore is not to be compared with the asymptotic, size independent, behaviour 
t~^/^ of the Potts model. The relevant point here is that u ^ u* > when 
t oo, if r > 4. Frozen states with density lower than in the closest packing 
have been observed recently also in different exclusion models (Grigera et al. 
1997, Fusco et al. 2001). 

Discussion and conclusions 

The results obtained for the exclusion model are very similar to those 
obtained for the Potts model. An important point to stress is that the 
two models are predicted to have the same symmetries and thus the same 
critical behaviour (Domany et al. 1977). However they are intrinsically 
different because there are features of the Potts model, hke surface tension 
and spontaneous magnetization, that are not shared by the hard particle 
model. Similarities in their behaviour support therefore the idea of = 4 
(and Tg = 4) as an upper bound for poly crystalline behaviour also in the 
Potts model. 

Ftom the above considerations it can be argued that these different be- 
haviours are actually related to the equilibrium phase diagram of the system. 

In fact, in two dimensions, for both Potts (Wu 1982) and hard particle models 
(Runnels 1983) the upper critical number of ground states is four. 

Another important point to remark is that the appropriate procedure to 
get the limiting energy is to take first the thermodynamic limit L oo, 
and then the infinite time limit t — > oo. This is a very important point 
that makes possible the observation of glassy phases in the kind of systems 
considered here. In fact systems with only short range interaction and no 
quenched disorder cannot have metastable states in the usual sense (Griffiths 
1964, Yeomans 1992). 

In conclusion we have shown that dynamics of simple lattice models ex- 
hibits a transition from (poly) crystalline to glassy behaviour when the num- 
ber of ground states increases. Investigations on the Potts model and on 
exclusion models support the hypothesis that the critical number of states 
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at which the transition takes place is related to the equilibrium properties 
and, for systems with this symmetry, is four. This work has been supported 
by grants from CCInt-USP (Brazil) and CNR (Italy). 
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Figure 1: Behaviour of the energy per site with time, for a g = 3 Potts 
model on a square lattice. Different curves correspond to different size (L) 
and boundary conditions {pbc = periodic, fbc = fixed). The dashed line is 
an extrapolation to t = cx) of the 1/ J{t) regime at large times. 
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Figure 3: The increase of the excess energy for increasing g, when t oo, is 
well visible from the time behaviour of the energy per site for fixed boundary 
conditionss (L = 1000). 
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